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SUMMARY OF RESEARCH 


Avalanche Diodes for the Generation of Coherent Radiation 

Work under this grant during the last period was concentrated m two 
areas: solid-state devices and characterization, and optimum imbedding 

I 

networks for realizing best performance. The BARITT device (Barrier 
Injection Transit Time Diode) has been under investigation, as reported 
below. In addition, other work previously reported has resulted in publi- 
cations, preliminary versions of which are attached as Electrodynamics 
Memos Nos. 28, 29, and 30. 

BARITT diodes are under investigation for possible application as 
microwave amplifiers and oscillators. Measurements have been made of 
diode noise figures m the frequency range of 4-6 GHz. Initial results 
indicate that a noise figure of 6-8 dB may be possible. Devices are under 
development for operation as amplifiers and sources in the frequency range 
2-4 GHz. 

We are investigating optimum device structure and fabrication tech- 
niques necessary for low noise performance. The possibility of cryogenic 
operation is being investigated, with a view toward improving the device 
noise figure. 

Models-for-thede vie efornonlinearand-small-signal-ope ration-are 

under investigation and optimal imbedding and deimbedding models are 
under consideration. 

A list of publications supported by NASA Grant NGL-22-009-337 is 
attached. 
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The lack of notation for derivatives in APL is an important 
barrier to its acceptance in many disciplines. This paper des- 
cribes a notation for derivatives and a procedure for accurately 
evaluating expressions containing derivatives. 

Operators _a nd F un ctions 

In describing APL syntax, it is frequently useful to distinguish 
functions from operators. Primitive scalar functions, primitive 
mixed functions, and defined functions are all "functions." The 
term "operator" is used for a svmbol which acts on one or more 
functions or variables, and produces, as a result, a new function. 

The new function then acts on its arguments and produces its 
results, an APL array. The operator, of course, must act before 
the function can, and therefore, the right- to-left rule for eval- 
uation in APL may have to be violated. 

There are three operators implemented in APL: reduction, 
inner product, and outer product. There has been suggested also 
a "scan" operator. These operators can, in principle, operate 
on defined functions as well as primitive functions, although 
they are not now so implemented. The derivative operator, to 
be discussed below, can likewise be used on both primitive 
and defined functions, but the implementation would probably 
be of little benefit if it were restricted to primitive functions. 


i^on s of _this__Not a 

There are three reasonable expectations for a notation for 
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derivatives. First, we might expect the notation to denote 
(and the implementation to calculate) the result of applying 
the new function produced by the operator, to an argument; 
that is, to denote the returned array. Second, if the purpose 
of such an operator is to return another function, one might 
wish the notation to denote (and the implementation to display) 
the resultant function. Third, we might expect numerical 
differentiation of functions defined only by a table of values. 

Of these three expectations, the second is difficult, both 
conceptuallv and practically , and the third is prone to numerical 
inaccuracy. The notation and implementation given in this 
paper fulfill the first expectation but not the others. 

One might also expect that it should be possible to devise 
a notation for integrals. The integral case appears to be 
much more difficult, however. 

Types__of _D e r iva_t i 1 ve s 

It may not be out of place to recall the wide varietv of 
derivatives that are encountered in mathematical ly based disciplines. 
These include the ordinary total derivative, partial derivatives, 
the substantive derivative of fluid mechanics, and the gradient, 
divergence, and curl in field theory, as well as high-order 
derivatives such as scalar and vector T.aplacians . Derivatives 
are taken with respect to any argument of a function, or with 
respect to any parameter. 

The proposed notation not onlv covers all these cases, but 
cats b •• further generalized. Differentiation is allowed not only 
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of functions v/ith respect to arguments, but also of functions 
or expressions with respect to global variables. 

Graphics 

Two now symbols are required for the derivative notation, 
one for an operator and the other for a function. The symbols 
fe and fe are used here, without anv implication that they are 
optimum. 

No ta t ion__f or Derivatives of Monadic functions 

Consider first a monadic primitive scalar function, 
with a scalar argument. The syntax for this function is Z+-F A; 
the derivative of F with respect to A, evaluated at the 
particular value of A , is denoted FM . For example, +b3 
is enual to l; is equal to "0.25; and L fe 3 . 5 is eoual to 

0 . 

Next, consider any monadic function, such as a primitive 
scalar function, a primitive mixed function, or a defined function. 
Again, the syntax is Z-*-F A. The derivative of F with respect 

_to A i s_denoted — FM-,— and— i-s— d e-f-i-ned— prov-id ed — A — and F—A— 

are conformable. The rule for conf ormabili tv is the same as 
the rule for the expression A+F A; that is, either A and F A 
have the same dimensions, or else one of them has only one 
element. If A is a scalar and FA is an array, then FM 
consists of the derivative of each element of F with respect 
to A, evaluated, of course, at the value of A given. For 
example, if the (velocity V of a particle and its position R 
are both functions of time T , and velocity and position in 
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3-space are represented by vectors of length 3, then the 
velocity is the derivative of the position, as expressed by 
the eouation ( V T) = RkT . 

If A and F A have the same dimension, then FM 
consists of the partial derivatives of each element of FA 
with respect to the corresponding element of A. In general, 
each element of the results depends on the values of all the 
elements in A t and is, of course, evaluated for the A 
in Question. This form of the derivative is useful, for example, 
in describing waveforms as a function of time: if FT is the 
function 2x* _ 5 *T which might be the natural response of a 
first-order physical system, then FVT would calculate an array 
consisting of the derivative of F evaluated at each of the 
times in the array T. 

Finally, if F is a scalar function of a vector • argument 
R which might stand for position in space, then FkR is a 
vector, each element of which is the partial derivative of F 
with respect to the corresponding element of R; that is FhR 
is the gradient of F. Thus, the relationship between electric 
f ie'l'd - vec"tor E and potential - PH'I ~ in electrostatics can be 
written in the form E - -PHIkR . As another example, 

( */ M ) = ( *//l ) *A except that it works properly even if OeA. 

In the last example, note that the reduction operator acts 
before the differentiation operator. 

Ou t e_r Der iv atives _o f Mo na dic F unctions 

It is freouently necessary to denote the partial derivatives 
of each element of the result of a function, with respect to 
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each element of its argument. This concept resembles that of 
the outer product, and a similar notation is suggested. Thus, 
for a monadic function F, F°.M is called the "outer deriv- 
ative". There are no conformability reouirements ; the outer 
derivative has dimension ( pF» . M 1 = ( pF A),pA. An example of 
the outer derivative is the velocity-gradient tensor in fluid 
mechanics, of which the rotation is the antisymmetric part and 
the strain rate is the symmetric part. Thus, if the velocity 

V is a function of space R, then the rotation is 

o . 5* ( Vo . b/? ) - . fe/? and the strain-rate tensor is 

0 . 5* ( V« . fe/? ) +$7° . fcF . As another example, + °.|?iF is the unit 

matrix of size N, N. 

Inn er Derivatives of Mon adic _Fu notions 

The inner derivative of a function F is denoted F+.M, 
and the conformability requirements are the same as for the 
inner product (F A)+.xA. As an example, one of Maxwell's 
equations in electromagnetic field theory states that the 
divergence of the magnetic field P is zero, or 0 =F+.feF. 

Generaliz ations of the inner derivative are possible by 

using scalar dyadic functions other than +. For example, 
F*.fc/I. For a function FT of time, its largest derivative 
at any time m a vector T is given by Ff.tT. 

I ' icfh -Ordxir. J2o r_i_v a_t_a_vo_s _of Mon ad i c Fun c t i on s 

Second and high-order derivatives of monadic functions~ 
are denoted similarly, with multiple use of the derivative 
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operator, for oxanple FbfeT. Thus, is the collection 

of all possible second-order partial derivatives of F t and 

has dimension equal to (p F A ) , ( p/1 ) , pA . 

An interesting second-order example is the Laplacian of 
a scalar function F of space coordinates R, which is Fo.fe+.feF. 
Mote that the outer derivative is performed before the inner 
derivative; that is, the derivative on the left before the 
derivative on the right, in contrast to the normal APL con- 
vention about order of execution of functions. If the function 
V R is a vector function of R f then the vector Laplacian 
is Vo .fe+.bF. As another example, consider Newton's law for a 
particle, where F is the force, M is the mass, and A is 
the acceleration, the second derivative of its position R with 
respect to the vector of times T : F-MxR^ feT. As another 

example, the homogeneous differential equation obeyed by a 
series RLC electrical network: 0 = ( LxJbfe2 , )+(FxJfe7’) + (J T)iC. 

Der ivativ es_of _Dy ad i_c JFun c t ions 

The notation for derivatives, inner derivatives, and outer 

_der.i_v_at;iv_es_of_dy_adic_func-tions_is_similar— to— that— for— monadic 

functions. Consider a dyadic function F with syntax Z*-A F B. 

The derivative of F with respect to B , denoted A FiB , 
must satisfy the conf ormability constraint suggested by the 
operation D+/1 F B . The result is the derivative (or partial de- 
rivatives) with respect to B or its elements, keeping A 
constant. The result is evaluated for the particular values 
of A and B presented. 
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A dyadic function can also be differentiated with respect 
to its left argument. In this case, the symbol t is moved 
to the left (that is, it appears between the function name 
and the independent variable) . The rules are similar and in 
particular, for the outer derivative, its dimension is 
( pA » A >F S) = ( pA) ,pA F B . 

As an example, if functions for current density J and 

charge density 0 are dyadic, with time on the left and space 

on the right, then the equation for conservation of charge, 

which appears in conventional notation as V*J + = 0 would 

0 t 

appear in the new notation as 0 = (T R. In a similar 

way, the substantive derivative, or "total derivative" of fluid 
mechanics ~ + V.V can be applied to a dyadic function N 

of space and time: ( Tkll R)+(T N» . R) + .*V. 

Some other interesting examples: The derivative of the 

maximum function A[ is equal to B^A . The derivative of 
the times function with respect to one argument is the other 
argument, that is (Ak*B)-B. If A is a vector then 
is a diagonal matrix with the elements of A along the main 
diagonal. 

Dc r i' va t i_ve_s _o f Kxpressjo n s 

The notation so far is suitable for differentiation of 
a function with respect to an explicit argument. It may be 
desired to differentiate a function with respect to a global 
variable, or to differentiate an APL expression with respect 

\ 
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to some variable in the expression. This can often be denoted 
by defining an extra function with the desired independent 
variable as an argument, but that technique is neither convenient 
nor natural. A slight generalization of the preceding notation 
is useful . 

The derivative of an APL expression E with respect to 
a variable X is denoted EkX . If the expression is more 
complicated than a single niladic function, it must be parenthesized. 
The variable on the right must be a variable name, perhaps 

indexed, but not a constant or an expression. The expression 
might involve the independent variable explicitly, as in the 
identity ( N*A *N-1 ) = ( A *N ) M ; or perhaps as a global variable in 
a function, as in the identity ( (AxB)kX)=((AkX)*B)+AxBbX; or 
perhaps in both roles . 

Outer and inner derivatives of expressions with respect 
to global variables are denoted similarly, E°.bX and E*-.^X. 

Note that since X must be a variable name rather than 
an expression, the following examples are syntax errors: 

EkX + 2-, Ek(X); Ek-X; EkX,0pC<-2; EbX+-3. 

Impl eme ntation 

The implementation of this notation appears to be 
straightforward, although perhaps demanding of space and time. 

The basic idea is that during execution of a function (or 
expression) being differentiated, a more complicated procedure 
is used. Every execution of a primitive function is accompanied 
by the calculation of the corresponding derivative. Thus, as 
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the defined function is executed, line by line, not only 
are the values of the expressions calculated, but also their 
derivatives, always with respect to the original value of 
the argument in question. When the function finally returns, 
the value of the return variable is discarded and only the 
derivative is used. 

This procedure, which is defined in detail below, in 
principle can work for any number of derivatives of arbitrary 
order . 

Details _of _t h e Implementation 

A new system-dependent function which I will denote here 
as i50, is required. Initially, this is set to the empty 
vector i 0 . 

Consider a function F whose argument A is the independent 
variable, with respect to which the differentiation is being 
carried out. When the execution of F begins, several things 
happen. First, the value for 150 ceases to be the empty 
vector and becomes the vector ,1. Second, the independent 
variable is accompanied by + that is, the outer derivative 

of /I with respect to itself. Other variables, both global 
and local, are also considered to be accompanied bv their partial 
derivatives with respect to A f in every case an array of 
zeros. Thus, all variables, both global and local, are considered 
to be not only their values but also their outer derivatives with 
respect to A . 

As the function is executed, line by line, each execution 
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of a primitive function is accompanied hv the evaluation of 
the derivative of the result with respect to A; thus, for 
everv array qenerated there is also qenerated another array 
consistinq of the outer derivative of that ouantity with respect 
to A. The chain rule of differentiation is used. Naturally, 
the derivatives so calculated will all be zero until A actually 
enters into the calculation. 

If another defined function is called within the function 
beinq differentiated, this mode of operation continues and when 
this function terminates execution, not only the value of the 
return variable but also its outer derivative with respect to 
A will be returned. 

A qiven differentiation is completed when the function (or 
expression) beinq differentiated finally returns, or else is 
removed from the state-indicator stack by the right arrow. When 
the function returns, the return variable has both a value and 
its outer derivative with respect to A. The value is discarded 
and the derivative is returned. In the case of the derivative 
or inner derivative, appropriate operations are performed to 
extract the desired derivative. At the same time, 150 is reset 
and all qlobal variables have their accompany inq derivatives 
erased . 

When more than one derivative is beinq considered, a similar 
process applies. High-order derivatives with respect to the 
same argument cause values of 2 or more to be inserted into 
ISO, rather than 1, and variables are augmented with both 
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the first and second derivative with respect to A. 

Derivatives with respect to other variables mav he 
encountered while anv given derivative is pendina. In that 
case, 150 will contain not only the 1 (or hiaher number) 
associated with the initial derivative, but also a number 
associated with the more recent derivative. Thus, 150 is a 
push-down stack, with one entry corresponding to each independent 
variable. During execution of functions, not only the derivatives 
with respect to each of the independent variables, but also cross 
partial derivatives have to be calculated. High-order derivatives 
of any order can be calculated recursively by using onlv the 
formulas for first-order derivatives. 

While a derivative is being calculated, that is, while 150 
is not emptv, the normal APL environment is changed in a few 
ways. This can be detected while the function is executing, or 
while it is suspended. First of all, the state of any pending 
derivatives can be interrogated by asking for i50; pi50 gives 
the total number of different independent variables under consid- 
eration, and +/I50 gives the highest-order derivative being 
calculated. Second, execution of the system commands )SI or 
)SIV will cause an extra symbol to be printed next to the 

names of functions being differentiated. Third, the derivative 
of any expression E with respect to the Nth independent variable 
is denoted by N^B , or for the outer derivative, A ,’o.fcF, or for 
the inner derivative, /<’+.££’. This feature allows functions 
to branch according to the values of derivatives, and also allows 
the display of derivatives for debugging. 
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The symbol b can also be used monad ica llv , as in 
or or + .fcF to denote the derivative with respect to 

the most recent independent variable. 

Needless to sav, much space and time will be saved by 
adopting the strategy of not actually storing derivatives that 
are identically zero. Thus, if execution is suspended while 
a derivative is pending, independent calculations can be 
carried out without necessarily calculating all the derivatives 
which end up eorual to zero anyway. 

It may sometimes be desired to obtain both the result of 
a given function and its derivative without re-executing it, for 
example in case the function uses terminal input or makes changes 
in global variables . This is easilv accommodated with the 
proposed notation. The last tine the return variable is assigned, 
also assign some alobal variable. When the function returns , 
the aloba] variable will retain the value of the function 
whereas the return variable will return the derivative. 

rnresolv ed_ P rob. lei m 

The derivative of every primitive function in APL can be 
written in terns of other primitive APL functions by use of 
the chain rule, with one exception. The exception is the scalar 
function . The derivative of the aamma function cannot be 
expressed in terns of other implemented APL primitives. Further- 
more , i t is not known whether derivatives of arbitrary orders 
can be calculated through a recursion formula. 
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Conclusions 

The notation for derivatives in this paper is convenient, 
yet general enough to cover probably all instances of differ- 
entiation in mathematics. It does not, however, perform or 
denote numerical differentiation, nor does it indicate the 
display of functions which are derivatives of other functions. 
The implementation outlined here is believed to be feasible, 
and if realized, would lead to a significantly more useful APL. 
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Introduction 


Several attributes of two-port networks are well known, 
including reciprocity, passivity, losslessness, unconditional 
stability, and symmetry. Not all two-port networks have all 
these attributes, of course, and two-ports can be classified 
according to whether each attribute is or is not obeyed. For 
example, we have the well known class of reciprocal networks, 
or unconditionally stable networks. Many of these attributes 
can also be defined for one-port or n-port networks . Each 
can be expressed as a constraint on the two-port parameters, 
for example reciprocity implies Z ^ = Z 21 and symmetry im- 
plies, in addition, Z^ « Z 22* Each of the above attributes 
can be associated with a numerical factor that is either in 
a certain range or equals a certain value if the attribute 
is present. Other information can sometimes be conveyed by 
the value of the factor even if the attribute is not obeyed. 

For example, the "reciprocity factor" Z 2l/ Z l2 e< ? ua ^- to ^ 
for reciprocal networks, and if it is different from 1, 
its magnitude suggests the degree to which the two-port is 
unilateral, with values of » or 0 for the unilateral case. 
Similarly, the "symmetry factor" Z n/ Z 22 ^ or rec ^P roca l two- 
ports is 1 for symmetric networks, but if it is not 1, its 
value indicates how the impedance level is transformed by the 
two-port. 



Tightness 


A new attribute, "tightness," is motivated by the obser- 
vation that the two networks in Figure 1 are somehow basically 
different, yet their difference is not suggested by any of 
the attributes listed above . (Note that they are both - 
reciprocal, passive, lossless, unconditionally stable, and 
symmetric.) 

A tight two-port network is, intuitively, defined here 
as one for which changes in the termination at the input or 
the output are greatly visible from the other port. A network 
that is not tight in this sense is one in which the output is 
more or less decoupled from the input, or vice versa. A use- 
ful numerical measure of the tightness of the network is found 
by comparing the input impedance Z^ n for two different load 
impedances. For the two loads it is convenient to choose 
open and short circuits, and the "tightness factor" TF is de- 
fined as the ratio of the open-circuit impedance to the differ- 
ence of the impedances. That is, 


TF 


oc 


oc 


- Z 


sc 


( 1 ) 


A simple evaluation gives this in terms of the Z or Y or ABCD . 
parameters as 



Parker, Peskin, and Chirlian^”^ have observed that the ratio 

Z /Z (and therefore the tightness factor) is independent 
oc sc 

of which port is regarded as the input. 

A tightly coupled network has a value of TF close to 1, 
and a network in which either the input or the output is de- ' 
coupled from the other has a tightness factor of <». This 
includes networks in which the input and output are uncon- 
nected, like Figure 1-b, as well as unilateral networks. A 
gyrator has a tightness factor of 0. 

The tightness attribute, like the others mentioned above, 
can be tested at any frequency, and the tightness factor can 
be evaluated as a function of frequency, either numerically 
or analytically. It is possible for a network to be tight 
at some frequencies , but not at others . 

Examples 

The tightness factors for several simple two-port net- 
works are given in Figure 2. Note that in the case of a mutual 
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inductor, the tightness factor is 1/k where k is the coupling 
coefficient. In the examples the tightness factors are all 
real, although in general complex values should be expected. 

7-10 

This concept resembles in some ways that of compactness 
but is slightly different, and can be applied to any linear 
two-port, whereas compactness is usually restricted to lumped, 
or lossless, or RC networks. To relate these two concepts, 
a network is compact if it is tight at all its natural frequen- 
cies . 

5 

Anderson and Ku have discussed -the ratio of open-circuit 
to short-circuit impedances for n-ports, and their approach 
might be useful as a basis on which to extend the tightness 
concept to n-ports . 



c «. j ~ ' ' a) 


Figure 1. Two extreme examples of two-port networks. 

Network (a) is tight whereas network (b) is not. 
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Figure 2. Values of the tightness factor for several 


two-port networks . 
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I . INTRODUCTION 

MARTHA is a notation for denoting electrical networks, 
and it is also a computer program which uses this notation. 

A. General Description 

In MARTHA, the notation, every network is either an ele- 
ment (resistor, capacitor, transmission line, etc.) or else 
one or two previously defined networks wired together. The 
basic idea is illustrated in Figure 1. In this example, S 
and P are "wiring functions" and R, L and C are "element 
functions". 5 and P are dyadic; P, L and C are monadic. 

The functions used in Figure 1 are sufficient to describe 
any series parallel network containing linear resistors, capa- 
citors, and inductors. For more complicated networks, other, 
similar, functions in MARTHA are used. 

As a notation, MARTHA is an alternative to the widely used 
schematic diagram. It can be used for communication from one 
person to another (for example, for documentation purposes), 
or from a person to a computer, or from a computer to a person. 
It is useful both in analysis programs' (where the user writes 
it and the computer reads it) and in synthesis programs (where 
the computer writes it and either the computer or the user 
reads it) . 

MARTHA, the computer program, is a network-analysis program 
imbedded in the interactive language APL , and the syntax 
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( C If - 6 ) P (i? 2 ) 5 ( L 0.01) 

Figure 1. Illustration of MARTHA notation for a model of a 
parallel tuned circuit. The functions R , L, and C define 
elements, and the functions 5 and P perform the wiring. The 
functions R, L, and C are "monadic" (having one argument, 
located on their right) , and 5 and P are "dyadic" (having two 
arguments, one on each side) . These five functions are suffi- 
cient to denote any linear series-parallel RLC network. 


- 5 - 


resembles that of apl in many v/ays . The simplicity and 
versatility of the input language make MARTHA relatively 
easy to use. It differs from most of the other programs 
in this book in that it is naturally oriented toward ports 
rather than toward nodes. In MARTHA node matrices are 
never calculated or used, so that time-consuming matrix 
inversion is avoided altogether. There are no difficulties 
associated with capacitor or inductor loops or tiesets, and 
it is never necessary to find the eigenvalues of any large 
matrix. Analysis is carried out at all frequencies simul- 
taneously, using APL's fast array handling. 

B. Capabilities 

The program MARTHA performs frequency-domain analysis 
of 1-port or 2 -port networks which are made up of 1-port 
or 2-port linear elements wired together so that at every 
stage in the construction, only 1-ports and 2-ports are used. 
Examples of networks of this type are most amplifiers, fil- 
ters, and microwave systems. Networks may be active or pas- 
sive, may be reciprocal or nonreciprocal, and may be lumped 
or distributed (or a combination of both) . 

Engineers may use MARTHA for circuit analysis without 
knowledge of APL. The simple and uniform notation in MARTHA 
is advantageous- to beginners, as is the interactive nature of 
MARTHA. . However, users who knov; APL can write their own func 
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tions to control MARTHA . For example, network definitions, 
including both parameter values and topology, can be varied 
under program control. Advanced users can write synthesis 
algorithms using MARTHA notation for the resulting network. 

In other words, MARTHA can be used as a programming language 
besides simply as a program. 

Although MARTHA is intended for general-purpose network 
analysis, different users can make use of rather specialized 
portions of MARTHA to make up their own special-purpose analysis 
and synthesis systems. 

Microwave engineers can use the distributed elements in 
MARTHA, and wave variables (with complex, frequency-dependent 
nornaliza tion if desired), and Smith-chart plottinc. Aspects 
of interest to amplifier designers include various transistor 
models, calculation of several measures of gain and stability, 
and MARTHA's ability to plot Nichols charts and the U/A "gain 
plane". Filter designers can make use of MARTHA's ability 
to scale frequency and impedance, and to perform high- pass and 
band-pass transformations. For active filters, there are 
four operational-amplifier models of differing complexity. 

Of interest to experimentalists is MARTHA 's anility to work 
with tables of measured performance, and use them for calcula- 
tion and interpretation as elements. Engineers concerned with 
model making will find an extensive repertoire of basic elements 
including 16 controlled sources, negative- impedance convertors, 
gyrators, and elements whose response goes with complex frequen- 
cy to an integral power. 
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The user can define many different networks, and wire 
the- o together or analyze them when desired. Two or more 
networks can be analyzed at the sane time, and the results 
compared. There are more than a hundred different resnonse 
functions that can be recuested. The results of an analysis 
can be printed, or plotted vs frequency, or vs any network 
parameter (on linear or log scale) , or vs another response. 

The response of one network can be plotted against the response 
of another network, or against a numerical table of values. 
Alternatively, the results can be stored for later display 
or calculation, possibly using other results. 

MARTHA incorporates an extensive set of tools for de- 
fining, editing, manipulating, and interpreting tables of 
numerical values. Such a numerical function of freouency 
( FOF ) can be interpreted as a quantity to be printed or plotted 
(alongside the results of a normal MARTHA analysis) or as the 
impedance, admittance, or scatterina coefficient of a numeri- 
cally defined element. < 

C ♦ Documentation 

MARTHA is described in the book, Paul Penfield Jr., 

"MARTHA User’s Manual”, The HIT Press, Cambridge, Hass., 

1971. That book completely covers versions of MARTHA 
dated 71. For versions of MARTHA dated 73 (and this in- 
cludes all versions on commercial time-sharing computers) 
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sevcral improvements arc described in the pamphlet, Paul 
Penfield Jr., "MARTHA User's Manual, 1973 Addendum, The 
HiT Press, Cambridge, Mass., 1973. Users of MARTHA can 
determine their version date by referring to the line which 
starts CIRCUIT ANALYSIS BY MARTHA at the top of each print 
or plot: either 71® or 73® appears in that line. 

Several other publications describing MARTHA or the ideas 
behind it appear in the bibliography at the end of this 
chapter. See references 1-4. 

Besides these publications, there is extensive on-line 
documentation in MARTHA; for information- on how to access 
this , type 

)LOAD 100 ROW MART HA 

DESCRIBE 

D. Availability 

MARTHA is available both from commercial time-sharing 
computer companies, and for use on separate machines that 
run APL. Inquiries should be directed to the Manager of 
Software Services, The MIT Press, 28 Carleton Street, Cam- 
bridge, Mass. 02142. 
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II. PROGRAM STRUCTURE 

In APL, programs are stored in "workspaces". MARTHA 
consists of nine workspaces of which one is purely documenta- 
tion, and seven constitute the "MARTHA library". The basic 
workspace* , 100 MARTHA, contains about 70 cooperatina "fore- 
ground" APL functions that the user may call directly, several 
"background" utility functions called by the "foreground" 
functions, and a few global variables. The foreground func- 
tions fall into six categories. First are functions which 
create elements, such as the functions R , L and C . Second 
are functions which wire networks together, for example S 
and P. Third are functions that calculate the response of 
the network, for example impedance, admittance, reflection 
coefficient, or VSVIR of one-port networks, or any of the two- 
port parameters or various gain or stability measures of two- 
port networks. Fourth are functions which can modify the 
response functions by taking the real part, imaginary part, 
magnitude, etc. Fifth are functions that help specify the 
format of the output. Finally, are some miscellaneous func- 
tions to aid in defining networks. 

The MARTHA library contains additional more specialized, 
functions in these same categories. It is a standard feature 
of the APL workspace-storage system that individual functions, 
or groups of functions, can be copied into the user's active 
workspace; by so doing the user can select those functions 
from the library that he needs and leave the rest behind, thereby 
creatine his own, personalized, version of MARTHA. 

* On some rcr outers, the number in the v,-orkspace nares is 
different frew 100. 
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ll I. NETWORK ELEMENTS 

Table I and Figure 2 show the network elements defined 
in MARTHA* The two-port networks WR and WTHRU are constants, 
and all the others are monadic functions . The simple functions 
R, L, and C (resistors, inductors, and capacitors) have been 
illustrated earlier. The function L, when used with an argu- 
ment of length 3 (that is, a vector with 3 numbers in it), 
produces a mutual inductor; thus the function L produces dif- 
ferent elements according to the length of its argument. The 
same is true of the functions TEM , WG, and OPAMP . The wave- 
guide function WG produces a length of waveguide if its argu- 
ment has 3 numbers in it, but if only 2 are present (the 
length is absent) the result is the frequency-dependent charac- 
teristic impedance of the waveguide. The vaveguide analysis is 
valid both above the cutoff frequency, where the characteristic 
impedance is real, and below the cutoff frequency, where it is 
imaginary. Examples: 

R 40 (resistor, 40 ohms) 

L .015 (15-mH inductor) 

L .015 .02 .01 (2-port mutual inductor) 

IT 3 (ideal transformer, turns ratio 3:1) 

WG 1 E9 377 (matched load of v/aveguide with cutoff frequency 
1 GHz, 377-ohm characteristic impedance at f = °) 

WG 1E9 377 .6 (60-cr length of vaveguide) 

WG 1E9 377 90 DEGREESAT 2E9 (quarter-wave section of guide) 

WG 1E9 377 .4 FOPDIEL 2.5 (dielectric-filled guide) 
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Table 

1. 

Elements defined in MARTHA . 

In the 

equations , 

5 is 

j 2 IT f . 

The 

library workspace 100 MARTHAE 

contains 

many 

additional 

elements . 



blbmert 

TIPS 

NAME 

ARGU ,JP NT VECTOR 

REOUIRED 

EQUA TIONS 

resistor 

1-PORT 

R 

RESISTANCE RES IN OP MS 


V=RFS=I 

CAPACITOR 

1-PORT 

/• 

CAPACITANCE CAP IN FAN ADS 


I=S=CAP=V 

INDUCTOR 

1-PORT 

L 

INDUCTANCE IND IN BENRIFS 

IND* 0 

V=S=1RD=I 

STRA IGRT-TRR0US3 
CONNECTION 

0- PORT 

WTHRU 

(NONE) 


V1-V2-. ll=-I2 

POLARITI REVERSE 

2-PORT 

VR 

(NONE) 


Vl = -V2i 11 = 12 

MUTUAl INDUCTOR 

2-PORT 

L 

INPUT S p LP -INDUCTA NC p L 1 IN HEN P IFS 
OUTPUT S p LF - INDU r TA NCE L 2 IN REN P IFS 

mutual induct AVCF m IN henries 

M*0 

V1=S={L1=I1 )+M*I2 
V2=S*(H*T1)+L2=T2 

IDEAL TRANSFORMER 

2- PORT 

IT 

TURNS RATIO N 

N*0 

V1 = R= V2 ; Il=-I 2*3 

operational 

amplifier 

2 -PORT 

OP AMP 

O p FN-CI ppr JIT VnirACv CATV A 
OUTPUT I* ,PP DA N r E ROUT IN OH"S 
INPUT IMPEDANCE RIN IN OEMS 

A* 0 
RIN* 0 

V1=RI3 *1 1 

V2 = (A=V1)*R0UT=I2 

operational 
a "P LITTER 

2-PORT 

O' 0 AMP 

0P p N • CZ prt JIT VOLTAGE CAIN A 
OUTPUT I t,p EDANrv p 0UT In 0 R**$ 

4*0 

1 1 = 0 

K2=M»l'l )+ROUT*I2 

0PER. AMPLIFIER 

2-POPT 

OP A MP 

VOLTAGE GAIN A 

A* 0 

11 = 0-. V2=A = V 1 

pield-eeeect 

TRANSISTOR MODEL. 

grounoed-sou=:e 

2-POPT 

P P T 

GATE-SOURCE CA P ACITANC P CCS IN FARADS 
GATE -DRAIN CAPACITANCE CGD IN FARADS 
TRANSCONDUCTANCE G '• IN ”HQS 

C u * 0 

n=5«(<r£;s«Ki )*ccd= vi-v2 

I2=(GM‘V1)*S*CCD*V2-V1 

BIPOLAR-TRANSISTOR 

model, crounded- 
EMITTrR 

2-POPT 

RTBRIDPI 

RESISTANCE RJ IN OH'*S 
RESISTANCE RPI IN Or MS 
CAPACITANCE C P T IN FARADS 
CAPACITANCE C"U IN FARADS 
TRANSCONDUCTANCE CM IN MHOS 

RPI* 0 
CM*0 

VI =VPI*RX=I1 

n=i vpi=(s*cpn*iRpn* 

S=CMU= VPJ- V2 

I2=(CM*VPI)+S*CMU*V2-VPI 

LOSSLESS TRANS- 
MISSION LINE 

2 -PORT 

TEH 

CHARACTERISTIC IMPEDANCE Z0 IN ORMS 
LENGTH LEN IN METERS 

20*0 

A==-J=02*LE3=P*3E8*DIEL • 
( V2-zo*i7)=A*vi*zo*n 
( V1-Z0=I1 )=A=V2+Z0=I2 

TRANSMISSION LINE 
CHARACTERISTIC 

impedance 

1-PORT 

TEH 

CHARACTERISTIC IMPED AN CF Z 0 IN ORMS 

20*0 

v=zo=i 

lossless WAVEGUIDE 
DOMINANT MODE 

2-PORT 

VC 

CUTO V F FR p OU p NCY PC IN HFRTZ 
IN FI NlTE-PREOU p NCY CPA RA CTE°ISTIC 
IMPEDANCE ZTN 9 IN OR' J S 
LENGT° LEN IN H P T PP S 

ZINP* 0 
~PCtP 

Z0=ZIRF*(1-(FC*P)=2)..S 
A=J=LE3=ZI3P=F*Z 0= 358 01' 
( V2-Z0 = I2)=( =-o2*A )«’,'! -Z 
(Vl-ZO*It)*( =-02= A )*V2+Z. 

UAVECUIDE CHARAC- 
TERISTIC IMPEDANCE 

1-PORT 

VC 

CUTOFF FREQUENCY PC IN HERTZ 

INF . -P pr Q. CHA p . I”P. ZIN p IN OHMS 

ZINP*0 
- FCtP 

V=I=ZIRP*( 1-(FC*F) = 2) = . 5 



VG FC , 71 t/F , LEU 

Figure 2. Elements defined in MARTHA . 


WG FC.ZT it 
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Koto the two auxiliary functions, DECPEESAT and FCRDIEL . 

The first is used to soecify an electrical lenoth in de- 
grees at a reference frequency (in this case 2 GHz) and 
the second is used for dielectrically loaded guides. 

The library workspace 100 PAR TEA E contains about 50 
additional elements and models. One of them is the voltage- 
controlled voltage source, VCVS. That is a 2-port element, 
with the controlling branch on the input and the controlled 
branch on the output. Current can be used as well as voltage 
for either the controlling or controlled branch, so there 
are three additional elements, named VCCS , CCV5 , and CCCS . 

In addition, if flux linkage or charge are possible as the 
controlling or controlled variable, there are 12 additional 
controlled sources, including for example a charge-controlled 
current source, OCCS . These are useful for modeling. 

Another function useful for modeling is ZPDE . If its 
argument is of length 2 , the first is an integer and the second 
is a coefficient, the result is a 1-port network whose impe- 
dance is equal to the coefficient times complex frequency s 
raised to the integer power. If the integer is 0, the result 
is a resistor; if it is 1, the result is an inductor. Values 
of the integer from -5 through 5 are possible. ZPDE also pro- 
duces, if its argument is of length 3, 4, or 5, a 2-port 
similarly defined element where the numbers in the argument 
are from the 2-port irpedance matrix. Other functions named 
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YFDE , HFDF , and A5CPFDE operate similarly. These are useful 
for modeling using a power series expansion. 

Other elements in the ’.M REE A library include gyrators, 
nullcrs, negative -impedance convertors, and attenuators and 
isolators, both for TEM lines and waveguides. Included also 
are functions for converting numerically defined functions 
°f freauency (FOF's) into elements. The FOF can be inter- 
preted as impedance, admittance, or reflection coefficient 
of a 1-port network, or as impedance, admittance, hybrid, 

ABCD, or scattering matrix of a 2-port network. 

The workspace 100 FA FT HAS contains, more elements than 
can be described here. Complete documentation appears in 
references 1 and 2. 

The library workspace 100 EARTH AX contains several 
auxiliary functions for working with MARTHA , many of which 
are useful models. Examples are functions to calculate 
characteristic impedance of coaxial or microstrip transmission 
lines, or calculate coaxial discontinuity capacitances, or 
calculate cutoff frequency and characteristic impedance of 
common waveguides. 
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IV. 'INPUT LANGUAGE 


I'ARTHA is interactive. The user sits at a computer 
terminal and types his input line, and gets an immediate 
reply. 

Some of the things a user must type have nothing to 
do with network analysis. After the user dials the telephone 
number of the computer which carries MARTHA (which may, of 
course, be many miles distant) he must log in by identifying 
himself. He must also load HART 'A A from the computer's public 
library, and perhaps copy seme of the I'ARTHA library and 
perhaps copy seme of his previous results, including pre- 
viously defined networks or models, or previously calculated 
results. At any time the user can save his work up to that 
point, or start over, or log off the computer. 

The rest of the user's time can be spent analyzing net- 
works. A network can be defined at any time, and an analysis 
can be requested at any time . The analysis might be of a 
previously defined network, or of a newly defined network , 
oossibly composed in part of previously defined networks. 

The general form of an analysis request in MARTHA is 


r pr i i:t 

PLOT 

J PLOG > 


SMITH • 
STORE j 


< output list> OF <network> 
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The first word used {PRINT t PLOT, PLOC , SPITE, or STOPS) fixes 
the basic foraat of the output; PLOG is used for a plot with 
a logarithrric scale for the independent variable. The form 
of the "ou tout list" is discussed in Section V. The word OF 
is required to separate the output list from the network 
description . 

The form of the network description is unique to MARTHA. 

Two types of APL functions are used to define networks . The 
first create elements, for example, resistors, capacitors, 
etc., and the second create networks out of other networks , 
by wiring then together. Element-definition functions were 
described in Scvuicn III. An example of a wiring function is 
S , which is dyadic and, like all dyadic APL functions, has an 
argument on each side. Thus, if A and B are 1-port (2-teminal) 
networks then A S B is the new 1-port network formed by putting 
the two networks A and B in series. Similarly, the function F 
wires two networks in parallel. Figure 1 shows a schematic 
diagram and a description in MARTHA notation of a simple series- 
parallel network . The functions S and P are sufficient to 
wire together all linear series-parallel networks. For more 
complicated topology, other wiring functions, designed to work 
on 2-port networks, are defined in MARTHA. 

There are fourteen wiring functions in MARTHA, including 
the functions s and P. These are shewn in Figure 3. To create 
2-port networks out of 1-port networks, the functions MS and P 
are used. These are monadic; that is, they have one argument 
instead of two . The argument appears on the right of the 
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function, as is true of all nonadic APL functions. Two 
2-port networks can be v;ired together several ways. The 
functions WPP, WPS, WSF, and WSS connect the two inputs 
and ths two outputs, using either parallel or series connec- 
tions at each port. These are useful in denoting feedback 
amplifiers; for example, the function WSS might be used for 
emitter degeneration. The cascade function WC is very com- 
mon. The ronadic functions WU and WROT are useful in convert 
ing grounded-enitter transistors to grounded-base or grounded 
collector configurations. Symmetrical filters can be denoted 
easily with WU. if the left-hand side of the filter is called 
A t then the overall filter is A' WC WU A. Finally, three tech 
nicues for converting a 2 -port network to a 1-port network 
by terminating the output port are shoim. The output port 
can be open-circuited (WTO) or short-circuited ( WTS ) or 
terminated in another network (wt ) • Note that WTO and WTS 
are monadic, but WT is dyadic, expecting a 2-port network 
as its left argument and a 1-port network as its right argu- 
ment. 

A precedence rule for the wiring functions must be 
established. By that is meant a rule for determining which 
of the functions are to be considered executed before others. 
For example, in ordinary algebraic notation, the expression 
AxB+3, written without parentheses, indicates that the multi- 
plication is to be performed before the addition. This is ar. 
example of the common rule that exponentiation is performed 
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bcfore multiplication and division, and those are performed 
before additions and subtractions. By way of contrast, the 
languace APL has a simpler precedence rule. The rule is than 
all functions have eaual precedence, and are executed strictly 
in the order indicated, from right to left, unless parentheses 
are used to delimit expressions which are to be evaluated 
first. Thus, in APL, Ay3+ 3 would be 4x(B + 3), rather than 
(4x3) + 3. The precedence rule for the wiring functions in MARTHA 
is the sane as that of APL. As a consequence, every monadic 
wiring function takes as its argument the entire expression 
to its right, and every dyadic function takes as its left 
argument the one object immediately to its left, and as a 
right argument the entire expression to the right. Paren- 
theses can, of course, be used in the usual way to surround 
expressions which are to be evaluated first. Parentheses 
are often required for left arugments of dyadic functions, 
but are never recuired (though they are permitted) for right 
arguments of functions . As an example, the expression 
Ai 3 A2 P 43 refers to the network of Figure 4(a), rather than 
the network of Figure 4 (b) , which would be written (A 1 5 42) PA 
Sere of the wiring functions expect 1-port networks and 
others 2-port networks as arguments. VJhat happens if the wrong 
kind of network is used? For example, what happens in A 5 5 
if B is a 2-port network? In order to allow all wiring function 
to oporate-on both 1-port and 2-port networks , two "auuoiratic 
conversion conventions" are adopted in VAPTHA. They are: 
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1. If a wiring function expects a 1-port network 
and encounters a 2-port network then the output 
is open-circuitcd and the input is used. This 

is eouivalent to using the function '/TO. 

2. If a wiring function expects a 2-port network 
and encounters a 1-port network, then the wiring 
function ,-P is automatically invoked to convert 
to a 2-port network. 

Several examples of networks are shown both in ' F TH A 
notation and in schematic diagrams , in Figure 5 . These in- 
clude the widely-used Darlington transistor connection, a 
half-ladder and a Vheatstone bridge. 

MARTHA incorporates several other functions that are not, 
strictly speakir.a, wiring functions, althouah they have as an 
argument a network and return as a result a network based upon 
the argument, All of these operate on both 1-port and 2-port 

networks, and return a network with the same number of ports. 

. 6-3 

The monadic function HAD converts a network to its adjoint , 
which is the network whose impedance and admittance matrices 
are the transposes of the corresponding matrices for the 
orininal network. The function ZSCALE is dyadic; its left 
argument is a number and the result is a network with the same 
topology as the net’-ork of its right argument, but with all 
elements scaled in impedance. Similarly, FSCALE performs a 
frequency scaling. These functions are useful in filter de- 
signs, where perhaps a 1-ohn, 1-Hz prototype is known. In 
a similar way, } I AVERT performs a low-pass to high-pass trans- 
formation , and F3P a low-pass to band-pass transformation. 
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Al S A 2 P A3 (A1 S A 2) P A3 


(a) (b) 


Figure 4 . Illustration of the precedence convention for 
wiring functions in MARTHA . 



(c) 


(d) 


Figure 5 (continued on next page) 
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100 Q 1 yF 
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51 


1 ] 

I 





Cl 


IN 


Q 2 


1 : 1 
U J~ 

;5 N 

r" 3 


-I 52 


(f) 


(g) 


Figure 5. Examples of HART HA notation. 

(a) Double-stub tuner with 50-ohm lines and adjustable line 
lengths LI and 52, and fixed separation 5 inches: 

C WTS TEH 50,51) WC (TEH 5 0 , 5 *. 0254 ) WC WTS TEH 50,52 

(b) Twin- tee filter, values in ohms and microfarads: 

((.75 R 2 0 0 ) WC ( C 2E~ 6 ) WC WS R 200 )WPP(WS C 1E~&)WC(R 100 ) WC VS C 15~6 

(c) Feedback amplifier: 

(Cl VSS R 500) WPP WS (P 1053) 5 (R 2053) P C 3E~6 

(d) Vheatstone bridge: 

((t75 R 40) VC R 50) WPS VR WC (WS R 20) WC R 30 

(e) Active allpass filter (reference 5) : 

( WS( R 100)5 C 15” 6 ) WPP( WS R 500 0 )WC ( ( 0PAHP ~1E6)WPP WS R 1000 )WC 
WS ( R 100 )P C 15"6 

(f) Darlington transistor connection: 

WR0T WROT (WR0T 02) WC WROT Cl 

(g) Half-lattice network: 

(WS 51) WPP (IT "1) WC WS 52 
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The latter tv:o ray be used in succession to dof ino a band- 
elimination filter. The function WCC returns a network with 
every element replaced by its "complex conjugate"; it is 
useful in defining ccr. j uratc-natch loads for optimum power 
transfer. Finally, the function '.'DUAL takes the dual of a 
network , changing topologv, and element values, and element 
types in the process . 

This last set of functions is useful v:hen I’ARTHA is 
used for design work . The y.APTHA library contains a function 
named UR AT 13 which orints network descriptions of any 
network. Thus , for example, if PPGTO is a 1-ohm, 1-JIz proto- 
type filter, then 

DESIGN 1C00 ZSCALE 156 F5CALF PROTO 
defines a 1000-ohm, l-::Hz filter. Then 

r 1 z; a rr — c* n r* C T r V 

will print its definition, and 

PLOT 03 IG OF DESIGN) , DB IG OF DESIGN WC WS R 800 
will analyze it (in one case after it is wired with another 
resistor) . 
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V. OUTPUT SPI’CIFICATIONS 


The general form of an analysis request in MARTHA is 


PRINT 

PLOT 

PLOC 

SMITH 

STORE 


< output list> OF < network > 


In this section the output list is described. 

The output list is, basically, a list of response functions 
of the network to be calculated. It also can include modifiers 
on the response functions, and some format requests. 


A. Hesponse Functions 


MARTHA can calculate over a hundred different response 
functions, of which the thirty in Table 2 are thought to be 
of wide interest, and the remainder of specialized interest. 

The response functions of Table 2 are in MARTHA , and the rest 
in the library workspace 100 MARTHAR. 

For 1-port networks, the response functions are the im- 
pedance Z, admittance I, normalized impedance Z and admittance 
l, reflection coefficient SC, and voltage standing-wave ratio 
VSWR. Of these, VSWR is real and the others are complex. If 
the network in question is a 2-port network, then (in accordance 
with the automatic conversion convention for wiring functions) , 
its output port is open-circuited and the input port is used. 
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Table 2. Response functions in MARTHA. Of these 30, 

2 6 are conplex and four real . Over 70 additional response 
functions are in the library workspace 100 MARTHAR . 



d/7 

MCA II' 1 G 

n RPR 71 r>5 O’’ 

n 

Lt 

C 

IMPEDANCE OF A 1-PORT HR T'.IQRY. 

prr ?*?f)RK 




V=Z*I 



Y 

c 

ADM ITT All CR OF A 1-PORT HE TWO PR 

NETWORK 




I- Ry-V 



SC 

c 

REF LECTI Oil COEFFICIEHT OF 1-POFT 

'IE TWO R Ji 9 

ZN 



B=SCxA 



Zll 

c 

I'lpRDA 11 CE ‘ V TR I X 

NETWORK 


Z-l 2 

a 

7 1 = ( Z 1 IxJl ) + (Z12 xJ2 ) 



Z21 

c 

72=(Z21xJl ) + (Z22 xJ2) 



Z22 

c 




Yll 

c 

ADM ITT A 'ICE 7, 'TRIX 

NET' 1 ORE 


Yl 2 

c 

11 = 1711x71 ) + f 712x1/2 ) 



Y21 

c 

12 = (72 lx 71 ) + (Y22 xI/2) 



Y2 2 

r 




till 

C 

HYBRID MATE IX 

NETWORK 


1112 

c 

71 = ( 711x11 )+(712x72) 



a 2i 

c 

12 = (921xJl) + ( 7/2 2x72 ) 

• 


1122 

a 




Sll 

c 

SCATTER IRC, ’L'TPTX 

*7 ■’? rn T v n v y 

it . ' .i . . 9 

ZFIN , Z’^7 

S 12 

c 

31= (Sllx/ll ) + (512 x/12 ) 



521 

c 

32 = (S2 lx/ll ) + (522x/2) 



S 22 

c 




ZI'l 

c 

IllPUT IMPEDANCE VI* 11 

UET rr OPK, 

ZL 

V T V 

c 

1 11 PUT AD 1 'ITT TWICE 11* 71 

TJ rj n r rp r> V ^ 

ZL 

SI •” 

c 

INPUT REFLECTION COEFFICIENT B1*A1 

HE TWO PE , 

ZL, Z N T 7 

ZOUT 

c 

OUTPUT 1'IPEDA UCE 

HE TWO R X , 

Z G 



V2*I2 VHF.il OUTPUT EXCITED 



TOUT 

c 

OUTPUT AD'II TTA11CE 

UE^l’O 9 

ZG 



1 2 v V 2 VP Ell OUTPUT EXCITED 



SOUT 

c 

OUTPUT REFLECT I Oil COEFFICIEHT 

NETWORK , 

ZG, Z "OUT 



B2*A2 VP Ell OUTPUT EXCITED 



VG 

c 

VOLTAGE GAI'l V2 *EG 

NET' 10 73 X , 

ZG, ZL 

AG 

R 

AVAILABLE GAIN POUT , AV * PI” , •* V 

11 r TWO oy i 9 

ZG 

IG 

R 

INSERT I Oil GAIN 

NETWORK , 

ZG, ZL 



pa U T( "F. TWO RX ) * PO U ? ( WTHRU ) 



PG 

R 

pomp GAIP -OUT* PI 7 ’ 

nr^’/n v i' 9 

ZL 

r* ^ 
i / 

R 

TP A US D 1J CF ? 0/ T 7 PO II T * PI ' 7 , A V 


ZG, ZL 



-26- 


In calculating SC, VSWR, Z, and Y a normalization impedance 
is required; in MARTHA the normalization irpedance for 1-port 
networks is known as Zll , and must be set by the user before 
analysis. It may be either an impedance (a real number) or 
a real or complex nurerical function of frequency ( FOF) or 
any 1-port network defined in MARTHA notation (in which case 
the impedance of the network will be used) . 

For 2-port networks, there is a wide variety of response 
functions. If a 2-port response function is requested of a 
1-port netv;ork, then (in accordance with the automatic conver- 
sion convention for wiring functions) the wiring function WP 
is assumed to be invoked. The response functions available 
include all the common 2-port parameters (impedance, admittance, 
hybrid, ABCD matrix, and others) and the corresponding matrices 
for wave variables (including the scattering matrix and scatter- 
ing transmission matrix) . For calculations of the vave-variable 
responses, normalization impedances at both the input and the 
output are recmirea; these are known as ZNIN and Z110UT, and may be 
different both from each other and from Z/7 . Each of these, like 
ZN , may be a real constant, real or complex numerically defined 
function, or any MARTHA 1-port network. 

Many of the 2-port response functions depend also on the 
generator and/or load impedance. A 2-port network is assumed 
to be terminated by generator and load as shown in Figure 6. 
the variables ZG , ZL , and EG may, like ZH , be specified in 
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several different rays, but they rust be specified before 
analysis . Examples of response functions include the input 
impedance, admittance, reflection coefficient, and VSWR 
(these all depend on ZL) and corresponding output Quantities 
(these all depend on ZG) . Other examples are various measures 
of gain, including power gain, voltage gain, open-circuit 
voltage gain, voltage ratio, transducer gain, insertion gain, 
insertion voltage gain, unilateral gain, conjugate-match gain, 
and available gain. Also available are various characteristic 
impedances of the netv’ork, includincr image impedance, itera- 
tion impedance, and conjugate-natch impedance, along with 
corresponding admittances, reflection coefficients, and pro- 
pagation constants. Also included are stability factors for 
amplifiers, and normalized impedances and admittances. Other 
response functions include the input and output voltages, 
currents, and wave variables. 

B. Modifiers 

The complex response functions appear normally in the 
form of real and imaginary parts. This may be changed by the 
use of modifiers. Table 3 lists the most important modifiers 
in MARTHA. For complex responses, the real part, imaginary 
part, magnitude, magnitude in decibels, angle, and phase delay 
are all useful. The appropriate modifiers are placed before 
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ZG 

WV 

EG INPUT 


GENERATOR 

Figure 6. Termination of two-port networks assumed when 
some of the response functions are calculated . The generator 
and load impedances ZG and ZL are not part of the netv. T ork 
definition. 



TWO -PORT 
NETWORK 


OUTPUT 


ZL 


LOAD 


Table 3. Modifiers in MARTHA . Of thses six, RE, IM , RAD, PD, 
and DEG are ignored if applied to a real response, and DB 
is 20*10 cMAGN ITUDE for complex responses and 10x10 <2ABS0LUTE VALUE 
for real responses. If no modifier is used, the real and 
imaginary parts of complex responses will result. Other 
modifiers are in the library workspace 100 MARTHAM . 


MODIFIER 

MEANING 

RE 

Feal Part 

IM 

Imaginary Part 

MAG 

Magnitude 

RAD 

Phase in radians 

DEG 

Phase in degrees 

DB 

Magnitude in dB. 

PD 

Phase delay 

REC 

Reciprocal 
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the responses in cuestion. Example: 

PR HIT HAG Z , DEG Z OF (R 1) PL .02 

Each of the modifiers acts on only the one response immediately 
to its right, contrary to the normal APL convention that func- 
tions have as their right argument the entire expression to 
their right. Thus parentheses are not ordinarily used in 
the output list. 

C. Formats 

Various format requests can also be inserted in the out- 
put lists, generally at any point in the list. For prints, 
the number of significant figures printed can be changed from 
its normal value of 5 by the function PLACES. Example: 

PRINT 7 PLACES DB IG OF CRY STA LFILTER 

For plots, several additional options are available. Plots 
are normally made against frequency as the independent variable, 
but if desired any one of the responses will act as the independent 
variable if it is preceded in the output list by VS. Alternatively 
if PAIRS appears in the output list, the first response function 
will be plotted against the second, and the third against the 
fourth, etc. Plots are normally 50 spaces wide and 50 lines high, 
to fit conveniently on one page (except those ma de by the function 
SMITH which are designed to fit a standard-size Smith chart) . 
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Other widths and heights can be specified in the output list. 
Plotting characters other than the standard ones can be 
specified by the function SYMBOLS. Normally MARTHA plots all 
dependent variables with different scales, selecting each so 
that the scales consist of round numbers, but still significant 
detail is shown in each plot. The dependent variables will 
all have the same scale (v/ith usually some loss of detail) if 
SS appears in the output list. The horizontal and vertical 
scales can be set to arbitrary values by the functions HSCALE 
and VSCALE. This is useful in magnifying certain critical 
regions of the plot; points falling outside the specified 
scales are simply ignored. 

These plotting format functions are illustrated by some 
examples : 


VLOG 'MD' SYMBOLS MAG Z, DEG Z OF NETWORK 
PLOT SS Zll, Z12, Z22 OF FILTER 

PLOT ~20 20 HSCALE "180 180 VSCALE DB RR VS DEG RR 
OF AMPLIFIER 

SMITH Sll, SIN, 512 OF FILTER 

PLOT PAIRS (Z OF AMPl), ZIN , Zll OF AMP2 

The third request prints the Nichols chart for the return ratio 
RR of the amplifier, and the fourth request produces a standard- 


size Smith chart 
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VI . MODELS 

A computation using MARTHA (or any other program) will 
only be accurate if reasonable care is used in modelling the 
network. Devising suitable models for a network is usually 
the most challenging part of any analysis. MARTHA cannot, 
of course, do the modelling job for the user, but it does offer 
him a selection of elements useful for models, the possibility 
of numerically defined elements , and the option of creating 
other elements in the form of "user-defined elements" . 

A. Tvoes of Models Allowed 

Models in MARTHA may be any network with topology definable 
in the MARTHA notation, using any of the built-in MARTHA ele- 
ments, along with any linear 1-port or 2-port element not already 
in MARTHA , for which the user is able to supply an APL algorithm 
for computing its impedance (for 1-port elements) or its ABCD 
matrix (for 2 -port elements) . 

Note that M.aRTHA does not distinguish models from elements, 
or from networks containing several elenents wired together. 

If any user linds that he uses a given configuration frequently, 
then he can write a simple APL function that returns that parti- 
cular network upon demand. For users who wish a model that is 
not representable by a network made up from built-in MARTHA 
elements, MARTHA allows user-defined elements without restriction 
as to complexity of the required algorithm. 
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Several elements of special interest in modelling were 
discussed in Section III. MARTHA can also handle numerically 
defined elements. To create one of these, the user types in 
a table of values and then calls one of the functions ZFOF, 
¥FOF, SFOF , HFOF , and ABCDFOF to interpret the table as numeri- 
cal values of impedance, admittance, or reflection coefficient 
for 1-port networks, or Z, Y, S, H, or ABCD matrices for 2-port 
networks. The resulting numerical models are then treated by 
MARTHA like other elements. Arbitrary frequency dependence 
is available this way. 

B. Built-In Model s 


MARTHA contains several built-in models, of which the 
most important are probably the hvbrid-ui and FET transistor 
models, and four simple models of operational amplifiers. 

These are considered as elements in MARTHA , and were dis- 
cussed in Section III. 

C. Innut Techniques and Format 

For use. r-aef ined models, a distinction should be made 
between those models that are networks of MARTHA elements, and 
those that are not. For the former, a simple APL function with 

an argument will usually suffice. As an example, consider the 

9 ------ 

Marcuvitz model for a window formed from a pair of semi- 
circular obstacles along the sides of a waveguide, Fiaure 7. 
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The eauivalent circuit shown has inductors with inductances 


L A = 


( ) 

2mf 2-rrfD 
c c 


L B = " 


16irf 


2 7Tf D . 
( — 2 -) 4 


Note that is negative . is the characteristic impedance 

of the waveguide at infinite frequency, f c is the cutoff fre- 
quency and c the speed of liqht. A function that returns this 

2-port network is 

7 P+VI'WOW A 

Cl] ;^300000000*( 6 . 28x/![ 1] x/l[ 3] ) 

[2] 5/W[2]x(,?*2)*(6.28x/ir 1] ) 

[3] L3+--*l 2]t( ( Y*4 )xl6x3. 14x/[l] ) 

[4] B-(W S L LB) UC ( L LA) l/C .75 L LB 
7 

where the argument A is a vector of length 3 containing the cutoff 
frequency, the impedance, and the diameter. A relatively small 
amount of APL programming ability is required to read or write 
this model. The normal APL function-definition scheme is used, 
and all APL primitive functions (and MARTHA functions) may be 
used . 

Next, consider the case where the model is so complicated 
that a network with MARTHA elements is not sufficient to describe 
it. The technique now is to use MARTHA'S user-defined -element 
capability. Two kinds of functions are required, the first to 
define the user-defined element, and the second to actually com- 
pute its response at each frequency. The first function is easv 
to write: it should return an APL vector of any length starting 
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Table 4. Required shape and interpretation of matrices 
returned by the user-written function EEEELErEET . 

Number of Columns 

Column 2 6 8 


1 

Re Z 

Re A 

Re A 

2 

In Z 

Im A 

Im A 

3 


Re B 

Re B 

4 


Im P 

Im B 

5 


Re C 

Pe C 

6 


Im C 

Im C 

7 



Re D 

8 



Im D 

Assumed 


(1+BC) 

D = 

A 


Properties 
of Result 

1-Port 

2-Port 

Reciprocal 

2 -Port 



Figure 7. Symmetrical window between two semicircular inductive 
posts in -a- waveguide and the rodel given by Marcuvitz' . __ 
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with the number 9. The rest of the numbers in the vector con- 
tain the parameters associated with the model. The second 
function is more difficult. it must be a monadic APL function' 
entitled NEUELEMENT . It is not called directly the the user, 
but instead will be called by one of the MARTHA functions when 
analysis is done. At that time the frequency is known (it 
need not be known during the definition of the network) and 
NEWELEMENT may refer to it. The function NEWEL EM ENT must 
return either a 2-dimensional matrix according to Table 4, or 
a MARTHA network, possibly incorporating such a matrix as an 
element. The first dimension of the matrix is in each case 
the number of frequencies, and the second dimension is 2, 6 or 
8, depending on whether the network is a 1-port network, a 
reciprocal 2-port network, or a nonreciprocal 2-port network. 

As an example, consider a length of coaxial transmission 
line with skin-effect loss. The TEM lines in MARTHA are loss- 
less, and the freauenev dependence of the skin-effect is not 
given exactly by any MARTHA element. For this example, the 

necessary parameters to specify the line are the inner and outer 
radii, and the length. A function that creates the user-defined 

element is very simple: 

V B*-CO?? r P L T ’’E A 
[11 3«- 9 2 . 5 2 7 7"7 

7 

The argument A is assumed to consist of the inner and outer radii 
and the length, all in peters. This function merely defines the 
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element; it does no: do any calculations. 

During HART HA analysis, the function NEWELEtiENT is called 
when this element is encountered. The real calculations are 
performed by that function. The ABCD matrix for this line is 
given by 

A = D = cosh HI 
B = Z o sin.h yl 
C = (sinh Y£) 

where Z is the length and 

Y * V( R + jojL) (G + jwC) 



and where the per-unit-length quantities L, C, R, and G are 
given by 3-0 

L 
C 

G 
R 


^ *n 
2tt r . 


2 irr. 


o 


in 


= 0 


= 2.5 2 x 10 " 7 /£ — ( — + — ) 


2tt r 


r . 

l 
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A function HE WELEKEHT that incorporates these formulas is 

V B .7 ~'UELr' f K! 1 T -1 ; CC\LL ; or> ;S 0 ;.Y0 ; ALBUM, ;BET^L 
C 1 ] rre-f-G . 2 8x3 . 3 5U7"12 v®4[4] T/1[ 3] 

[ 2] LL+- 2P 7 xe.'[ 4 ] v ■;[ 3] 

C 3 ] HR*-,' [2]x(,r*0. r j)x((T i 4[4]) + T/l[3])*6.23 

r 4] S0^(( (LL*CC)+( ( ( .05^0+2 ) + ( CCx 6 . 2 Ox ,F)*2 )■*<). 5)*2)*0. r > 

[5] ZO«--SB :-2x/?Ox(7Cx5 . 28x,F 

[6] /! L ? 7/: Z> «-/1 [ 5 ] x ? 7 v 2 x p 0 

[ 7 ] BETA.L*-/' [ 5 ] x 7 0 x SO 6 . 2 8 x , F 

[8] fl+((p,P), G)pO 

[9] 3[ ; 1 2 ]«-{?( 2 1 ° . oBFTAL ) x 6 5 o.oyt£P7/l£ 

CIO] 5[ ; 3 4 ] ( 2 1 » . oBFTAL ) x 5 6 °.o/!LP/7/.L 

Cll] PC; 5 6 ] -<- (SC ; 3 4]xpo,Cl.5] SO)+S[; 4 3] xZO , C 1 . 5 ] - XO 

C 1 2 ] SC ;5]«-?C ; 5 ] * ( Z 0 * 2 ) +S 0 * 2 
C 1 3 ] SC ;6]«-SC ;6]*(Z0*2)+S0*2 

C 1 4 ] SC; 3 4XPC; 3 4]x/?0,Cl.5] SO ) - SC ; 4 3]xZ0,Cl.5]-Y0 
7 

This function is not trivial, and some familiarity with fi.PL 
is necessary either to read or write it. 

The number 2.52£ - 7 in the function COPPERLIKE is a 
measure of the resistivity of copper, and is used in line C3] 
of HEWELEHEtlT . By including it as part of the element definition, 
other element-definition functions can use the sane function 
tJ EWE LEM E 11 T for other materials, for example 


7 B+mASSLIUC A 
Cl] 5«- 0 5. 017“ 7 , A 
7 


An important aspect of every model is its range of validity. 
This includes permissible values for both parameters associated 
with the model, as well as allowed freguency ranges. Users are 
advised to have their model -mat ing functions check the parameter 
values to see that they are reasonable. To check on frequency 
range, l .'ART BA can be made to print a warning whenever one or more 
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freauencies is outside the permissible range. To set the 

permissible range, the network is operated on by the function 

FLIUITS, which is a dyadic APL function with left argument a 

pair of frequency values limiting the allowed range. For 

example, the model for waveguide window given above is only 

valid (according to y.arcuvitz ) between f and 3f . If the 

c c 

function WINDOW had the following fifth line 

[5] B + <A :i], 3 x/lCl DFLIMITS B 

then whenever that network is analyzed, if any of the fre- 
quencies is outside that range, that fact will be reported, 
but analysis will continue. 
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VII. EXAMPLES 

Three fully worked examples are given in this section. 

The first is a transistor amplifier which is also treated in 
other chapters of this book. The second is a crystal filter, 
and the third a coaxial low-pass filter. The lines typed by 
the user are identified by being indented six spaces. The 
computer response generally is not indented. 

A. Transistor Amplifier 

Figure 8 shows the circuit diagram for the amplifier, 
together with the model to be used for the transistor. The 
bias source has been omitted. Eecause of the versatility of 
the technique of defining networks in MARTHA, the am.plifier 
may be defined in several ways. One way which is probably 
appropriate for this network is suggested by Figure 9, where 
the input and output coupling networks, and the emitter cir- 
cuit, are shown as seoarate two-port networks. Those circuits 
are straightforward, and are defined first (see EXAMPLE 1). 

The hybrid-pi transistor model used in MARTHA is generally re- 
garded as the best linear model at high freauencies. This 
model (see Figure 2) can be made to coincide with the transis- 
tor model to be used, if PX , CPI, and CMU are all set to zero, 
and CM is DETAiRPI. Then the output resistor R 5 must be added 
separately. The transistor is defined using this approach, and 



-40- 



Base Collector 



Emitter (fc) 


R 1 

50 ohms 

B 

98 


R 2 

20,000 ohms 

Cl 

10 .7 

pF 

R 3 

325 ohms 

C 2 

10 

yF 

R 4 

25 ohms 

LI 

1.7 

mH 


R5 2 Megohms 

(c) 


Figure 8. (a) Transistor amplifier. (b) model for transistor, 

(c) component values. 



INPUT ENCKT OUTPUT 


Figure 9. Transistor amplifier shown as four subnetworks 
wired together. 
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note that it is not necessary to separately calculate BETA +RPI , 
since the formula can be typed in and the calculated result 
used as part of the argument for the function HYBRIDPI . 

Refer now to EXAMPLE 1. First the workspace containing 
l-ARTRA is requested fron the APL public library. Then, the 
four subnetworks are defined, and then wired together and 
given the name x.vpi . This completes the network definition. 

Before analysis, the freouency vector F must be specified. 
A study of Figure 8 reveals no mechanism for limiting high- 
freauency response (this is perhaps an unrealistic circuit) . 

To view as wide a frequency range as is necessary, a logarithmic 
frequency sweep is used. The frequency vector F is set to 10 
raised to the power (the asterisk is used in APL for power) of 
a vector ranging frcn 0 to 11; this covers all interesting 
ranges . 

Next, a choice of resnonse function must be made. In 
this example, the voltage gain VC is most logical, and the mag- 
nitude (in dB) and phase (in degrees) are requested. The out- 
put recuest begins with the word PR I FT to set the basic format, 
and this is followed by the output list, then the word OF and 
finally the network definition. As expected, the voltage gain 
saturates at high frequencies at the rather high value of 138 dB, 
and the phase approaches -180 degrees. 

These results rust be taken with a grain of salt. Like 
all other network-analysis programs, ."ART 1 : A can give reasonable 
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answers only if reasonable models are used, and in the present 
case the circuit is probably not well modelled above 1 MHz. 
This fact does not destroy the usefulness of the circuit as 
an example of how to use MARTHA , but it does imply the need 
for better models . 

As an illustration of the use of MARTHA in judging 
different models, some rather simple improvements are made. 
First, a stray capacitance is placed across the inductor LI, 
and the resulting network named AMP 2 . This is, of course, 
still a crude model. Second, the transistor is re-defined 
with what might be typical values of RX , CMU , and CPI. The 
new network is called AMP3 . In each case, the network defini- 
tion is easy because the previously defined subnetworks are 
used. 

The voltage cains in the three cases are compared by 
simultaneously analyzing and plotting the gains of the three 
networks. To cover such a wide frequency range, a logarithmic 
plot is necessary, so the function PLOG is used (the function 
PLOT produces a linear plot) . Mote that MARTHA normally auto- 
matically selects the scales, but is here requested to use 
the sane scale for the three different responses. 
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B. Crystal rjlt^r 


This crystal filter, Figure 10, is adapted from a filter 
shown to me privately by Mr. William B. Lurie, who ascribed 
the design to Prof. G. Szentirmai. There are two stages, each 
with a half-lattice construction which is modelled with the 
aid of an ideal transformer. The four crystals all have the 
sane equivalent circuit, with different element values. The 
filter has a passband about 4000 Hz wide, between 8 MHz and 
8.004 MHz. It is designed to operate between source and load 
of 500 ohms. PvOfer to EXAMPLE 2. First M. ART HA is loaded, 
and then the four crystals are defined, next the two stages 
and the overall filter are defined. Next, the generator and 
load impedances ZG and ZL are set to 500 ohms, and the fre- 
cruency vector set so as to encompass the passband. In speci- 
fying the frequency vector, the index generator i was used. 

A plot 30 lines high of the insertion gain IG expressed in dB, 
and the phase of the voltage gain in degrees is reauested. 
Next, to see the ripple in the passband a little more clearly, 
the horizontal scale is deliberately set to the range from -2 
to 0 dB, with the aid of the function HSCALE from the MARTHA 
library. This request overrides the normal MARTHA practice 
of automatically setting the scale. Note that when a network 
is to be analyzed a second time for the same freouency vector, 
the previous .results, saved under_the name SAME, can_be used. 
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C 2 



Figure 10. (a) Crystal filter. (b) model for the crystals 

XI, X7 , X2 , and Z4. (c) component values. The element 

values in the crystal models are not listed since they 
appear in EXAMPLE 7. 
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Next, it is desired to inspect the transition region 
and the stopband, so the frecuency vector is redefined to 
cover the upper half of the passband and 4000 Hz above. 

The insertion gain, phase, and the magnitude of the input 
iir.pedance, are printed. 

An interesting phenorenon is observed at approximately 
8.007 nil z . The insertion gain passes through a peak and the 
phase changes by 180 degrees. Apparently there is a trans- 
mission zero somewhere in the vicinity, and to view this in 
more detail, the frequency vector is redefined so as to expand 
the range between 8.0068 MHz and 8.0070 MHz. The depth of 
the notch in insertion gain (105 dB) and still rather sharp 
transition in phase suggest that the zero is very close to the 
ju axis, and it is interesting to expand the passband again, 
between 8.00691 MHz and 8.00692 MHz. Note that the printing 
will not resolve these small changes in frequency so a larger 
number of significant figures is requested with the aid of the 
function PLACES from the MARTHA, library. The frequency sweep 
is now in 1-Kz steps, and the high resolution available in 
this analysis is due to the fact that MARTHA (like APL) uses 
double-precision arithmetic. 
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C. Coaxial Lc/7-Pnss Filter 

The third example, Figure 11, illustrates MARTHA'S ability 
to work with distributed as veil as lumped networks . The original 
design of this seven-section low-pass filter v»as done by R. Levy 
and T. E. Rozzi^, who stated that this filter is not particularly 
good except as an example of their design method, which they then 
used on a more practical, 23-section filter. Only the seven- 
section filter is analyzed here. 

The coaxial discontinuity capacitances of the filter are 
important to its design. These capacitances are models which 
should be used at every junction between conductors of different 
size, as shown in Figure 12. MARTHA has a function named 
C0AXD1SCAP in the MARTHA library which calculates the values 
of these capacitors, when supplied with the appropriate radii 
of the line. Like all lengths in MARTHA, the radii must be 
given in meters, and so a conversion from the dimensions in 
inches in Figure 11 is necessary. In EXAMPLE 3, the first four 
discontinuity capacitors Cl through C 4 are defined (the other 
four are equal because the filter is symmetric) , and the con- 
version from inches to meters, and from diameters to radii, are 
done during the definitions (multiplication by 0 .0254 meters 
per inch, and 0.5) . 

Next, the individual lengths of transmission line LI 
through Lb are defined. In MARTHA a length of line is specified 
by its characteristic impedance and physical length. In our 



Figure 12. Discontinuity-capacitor model for the junction 
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case we need to calculate the characteristic irrpedance from 
the dimensions. The function COAX in the MARTHA library is 
copied, and then used to do this . It expects as an argument 
a vector consisting of the inner and outer radii, and length, 
all in meters . 

Since the filter is symmetric, it is possible to define 
the left side (consisting of the first three lines with the 
discontinuity capacitors placed at all junctions) , and then 
use this definition twice in the final definition of the filter 
Making the definition in this way eliminates the need to enter 
the same numerical data twice, and therefore tends to reduce 
keyboard errors . 

The filter is designed to operate from a 50-ohm generator 
and into a 50-ohn load, so the generator and load impedances 
ZG and ZL are set to 50 . The VSUR at the input is wanted, and 
this response function VSWRIH is copied from the MARTHA library 
The freouency vector is set, and a plot is reouested. The 
result, of course, in in agreement with the analysis by Levy 
and Rozzi' 1 '^. 

Filters of this sort are periodic, and have other pass- 
bands. In this case there is another passband centered around 
12 GHz, so the stop band extends from 3 to 9 GHz. Levy and 
Rozzi stated that this seven-section filter had poor stopband 
attenuation, and the plot of the insertion gain up to 10 GHz 
reveals just how true this is. The maximum attenuation is 
only about 4 dB. 
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VIII. LIMITATIONS 

Every net' - or!:-ann lysis prograr has two types of limita- 
tions not discussed so far. One has to do with computer re- 
source limits, and r he other with problems that are ill con- 
ditioned . 

A. Resource Limitations 

Many programs have limitations on the number of nodes, 
or number of elements, or number of network defintiicns allowed. 
These limitations arise from the amount of space set aside for 
certain variables. MARTHA has, for practical purposes, no 
such precise liritaricns since it is imbedded in APL , which 
uses dynamic storage for all variables. The resource limita- 
tions of MARTHA are somewhat more difficult to describe. 

The primarv limitation is caused bv the finite size of 
APL workspaces (usually 32K bytes, each byte being 8 bits) . 

The MARTHA functions require about 1SK, leaving about 13K for 
the user's network definitions and temporary storage of re- 
sults. If possible, MARTHA, automatically analyzes with all 
frequencies at once; if this is not possible, it automatically 
selects a smaller number of frequencies, and repeats the 
analysis as many tires as necessary, until all freouencies 
are accounted for. However, even this can be insufficient in 
some cases with lengthy output lists. 


Experience has Shown 
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no trouble for analysis with up to fifty frequencies, provided 
not too r.uch of the available space is taken by unrelated 
functions and data. The APL error message US FULL indicates 
there is a problem of this sort. Sone of the computers that 
Carry l-'.APTF.A have larger APL workspaces , in v.’hich case this 
ceases to be a real problem. If it is necessary to analyze 
with a large number of frequencies and available space cannot 
be created by erasing unnecessary objects, then HARTHA has a 
function entitled ATATI'lE which can be used to repeat the ana- 
lysis for a small number of frequencies. 

B. Ill-Ccr.di tioned networks 

Every network-ana lysis scheme has its own set of networks 

for which it performs poorly, or perhaps not at all. During 

analysis, FAR TEA represents 1-port networks by their impedances 

and 2-port networks by their ABCD matrices. Thus MAPTHA is 

unable to handle networks for which these representations do 

not exist. For 1-port networks , this is only open circuits 

(or because of overflow or underflow problems, networks with 

37 

impedance magnitudes greater than about 10 or less than 
- 37 

10 ohms) . For 2-port networks, this is any network for 
which the output voltage and current are related by one equation 
not containing the input voltage or current. This class of 
networks are those_which produce no effect at the output when 
the input is excited, that is either those with disconnected 
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inputs and outputs, or those that are backwards unilateral. 
Generally \;hen such a situation exists, it is easy to rede- 
fine the network so as to avoid the ill-conditioned case. 

Certain elerent values, particularly zero, can lead to 
ill-conditioned networks, and Cable 1 lists several such cases. 
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IX . EPROR DIAGNOSTICS 

MARTHA recognizes tv7o classes of errors. One is errors 
so serious that MARTHA cannot reasonably proceed, and in the 
°ther case only sene of the calculations nay be in error and 
MARTHA continues . 

For non-fatal errors, MARTHA prints a warning but continues. 
Examples include use of freauencies outside the range specified 
by the function FLIMITS and frequencies at which by coincidence 
the network is ill-conditioned. Generally in the latter case 
trouble is encountered when a denoninator vanishes? the vrarning 
message is ATTEMPT TO DIVIDE BY ZERO. This also is printed when 
a user asks for a response function that is actually infinite, 
such as the admittance of a 0-ohn resistor. After this ressage, 
some of the results, but. not all, may be in error. 

Fatal errors are caught either by MARTHA or by APL . 

Among those caught by MARTHA are the wrong length for an argu- 
ment for an element-definition function. Errors caught by 
MARTHA induce an explanatory messarre, usually indicating what 
was expected, followed by the ressage MARTHA ERROR and a diag- 
nostic arrow pointing to the place where MARTHA got into trouble . 
These errors may often be regarded as a useful prompt, because 
by deliberately cornitting them, a user can remind himself 
what is expected. A complete list of these error messages 
appears in reference, 2. _ __ 
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Errors not caught by MARTHA, but rather by APL , are 
reported by the ressace RARE ERROR where NAME is the nare 
of one of the functions in MARTHA. This message may or may 
not be informative , but it is followed by an indication of 
the place v;r.ere MARTHA cot into trouble. 

To recover frcn an error a user should type a single 
right arrow -> followed by a carriage return, and then proceed 
after correcting the source of the error. 


-63- 


X. REFERENCES 

This chanter has not covered all aspects of using 
MARTHA ; in particular, no attempt has been made to describe 
the MARTHA library fully, because it periodically gets 
additions. Complete instructions for using MARTHA appear 
in references 1 and 2 belov/. Each workspace in the MARTHA 
• library has a variable entitled DESCRIBE which crives an 
up-to-date description of the contents. A succinct summary 
of MARTHA usage is available on-line in the vorksnnce 
100 H GUI' ART HA . 
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